CI Cam: A Shell-shocked X-ray Nova 
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ABSTRACT 

We present radio imaging observations of the 1998 outburst of the peculiar emission line star 
CI Cam, taken ~1, 4, 75, 82, 93, 163, and 306 days after the beginning of the 31.64 March 1998 
X-ray flare. The first two epochs show a resolved but compact (no larger than 12 milliarcseconds) 
radio source which becomes optically thin at frequencies higher than 5 GHz. The spectrum and 
brightness temperatures are consistent with synchrotron self-absorption, although free-free ab- 
sorption may also play a role. The later images show a large (120-350 milliarcseconds) oval-shaped 
or double ring remnant. The radio spectrum combined with the high brightness temperature in- 
dicates that the emission is synchrotron, while the morphology suggests that this is powered by 
a decelerating shock moving through dense circumstellar material produced by a strong stellar 
wind. The radio images of CI Cam are equally well fit by an expanding ellipsoid or two expand- 
ing rings; the former gives w 4.2(i — 50904. 1) 77 , with the major axis in mas and t the 
Modified Julian Date (MJD). The corresponding expansion speed in the plane of the sky was 
~ 12, 000 km/s over the first few days (for an assumed distance of 5kpc), slowing by a factor ~ 3 
by the time of the last observation almost a year later. The radio emission from all other X-ray 
binary transients has either been unresolved, or taken the form of highly collimated relativistic 
jets. We suggest that CI Cam represents a rare case where these jets were smothered early on by 
the unusually dense circumstellar medium. In this model CI Cam is the analogue to extragalactic 
supernovae formed by the collapsar mechanism, while the more usual X-ray binaries with rela- 
tivistic jets are analogous to the jets which escape those supernovae to form a subset of 7-ray 
bursts. 



Subject headings: binaries: close — stars: individual: CI Cam 
X-rays: binaries — X-rays: bursts 



radio continuum: stars — stars: flare 



1. Introduction 

XTE J0421+560 was discovered by the All-Sky 
Monitor (ASM) on the Rossi X-ray Timing Ex- 
plorer (RXTE) on 31.64 March 1998 (Smith et al. 
1998) as the first bright X-ray transient in several 
years, reaching almost 2 Crabs within a day of 
the first detection. Corresponding radio emission 
was detected soon after (1.9 April; Hjellming & 
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Mioduszewski 1998). Optical observations (Wag- 
ner & Starrfield 1998) identified the radio and 
X-ray transient as the star CI Cam, based on po- 
sitional coincidence with the radio emission, and a 
concurrent optical outburst. It was also detected 
by the Burst and Transient Source Experiment in 
the low-energy Gamma-rays, reaching a flux of 
1.6mCrab in the 20-40 keV band (Harmon et al. 
2004). 

CI Cam is a peculiar bright emission-line star 
with unusually strong Fe II lines (Downes 1984), 
indicating a very hot star or stellar shell. It also 
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MJD-50903.64 

Fig. 1. — The radio light curve from the Green 
Bank Interferometer (GBI), as produced by R.M. 
Hjcllming (priv. comm.). The green triangles and 
red circles represent 8.3 GHz and 2.25 GHz data 
respectively. After day 10 the data are averaged 
over 1 day and the error bars are the minimum 
and maximum for that day. The arrows show the 
epochs of our 7 imaging observations. 

has a near-IR excess (Dong & Hu 1991; Allen 
& Swings 1976). CI Cam has a rich and unusual 
optical spectrum and has been identified as a sym- 
biotic star (Chkhikvadze 1970), a Herbig Be star 
(The, de Winter, & Perez 1994), a B[e] star (Bel- 
loni et al. 1999), and a supergiant B[e] (sgB[e]) 
star (Robinson, Ivans, & Welsh 2002; Hynes et 
al. 2002). The recent consensus seems to be 
that it is either a B[e] or sgB[e] star. If it is an 
sgB[c] star it would be the only one identified in 
our Galaxy, all others having been found in the 
Magellanic Clouds. Its distance has been esti- 
mated to be anything from 1 kpc (Chkhikvadze 
1970) to 8 kpc (e.g., Robinson, Ivans, & Welsh 
2002) based on luminosity, spectrum, interstel- 
lar absorption, and kinematics. Here we take 
5 kpc as a reasonable compromise. CI Cam was 
not previously a known radio source, with two 
10.7 GHz observations in 1973 giving upper lim- 
its of 5 mJy (Altenhoff, Braes, Olnon, & Wendker 
1976; Woodsworth & Hughes 1977). 

CI Cam's X-ray /radio/optical flare puts it into 
another class, that of X-ray novae. If CI Cam is 
at a distance greater than 2 kpc then its X-ray 
luminosity indicates that there is probably a com- 
pact object (black hole or neutron star) in the 



system (e.g., Belloni et al. 1999). Orlandini et 
al. (2000) suggests that the outburst was the re- 
sult of thermonuclear runaway on the surface of 
a white dwarf, which would require a closer dis- 
tance. Even with all the different stellar identifica- 
tions for the optical source, it is clear that CI Cam 
must be a high mass X-ray binary (HMXB). How- 
ever, CI Cam's X-ray behavior is unique among X- 
ray binaries. For instance, unlike other HMXBs, 
CI Cam does not have persistent highly variable X- 
ray emission, and has a soft X-ray excess in qui- 
escence (Robinson, Ivans, & Welsh 2002). The 
persistent soft X-ray emission can be attributed 
to the sgB[e] star and the heavily absorbed harder 
emission may be an obscured compact object, but 
this is not certain, since this type of emission is 
extremely unusual for X-ray transients (Robinson, 
Ivans, & Welsh 2002). The flare itself peaked first 
in the X-ray, then in the optical, and finally at ra- 
dio wavelengths (Frontera et al. 1998), with sev- 
eral short, soft X-ray flares within ten days of the 
outburst (Frontera et al. 1998; Ueda et al. 1998), 
and rapid changes in the X-ray absorption column 
over the first few days (Belloni et al. 1999), as well 
as years afterwards (Parmar et al. 2000). The ra- 
dio behavior is fairly typical of X-ray transients, 
with an optically thick rise followed by an opti- 
cally thin decay (Clark et al. 2000), but with no 
sign of further flaring activity. An initial report of 
a relativistic radio jet (Hjellming & Mioduszewski 
1998) has since been shown to be an artifact due to 
a bad calibrator (Rupen, Mioduszewski, & Hjellm- 
ing 2003). 

In this paper we present a radio imaging study 
of CI Cam from a few days to roughly a year after 
the X-ray outburst, using the National Radio As- 
tronomy Observatory's (NRAO) Very Long Base- 
line Array (VLBA) and the United Kingdom na- 
tional facility Multi-Element Radio Linked Inter- 
ferometer Network (MERLIN). 

2. Observations 

CI Cam was observed on 1 April, 5 April, 
14 June, 21 June, 2 July and 10 September 1998 
with the full 10-element VLBA with the addition 
of one antenna of the Very Large Array (VLA), 
except for the 5 April 1998 observation, which 
was missing the Pie Town, Kitt Peak, and Owens 
Valley VLBA antennas. All of these observations 
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Table 1 
Observations 



Epoch 


days since detection 






frcq. 


time on CI Cam 


int. flux 


rms 


M J D a 


date 


of X-ray flarc b 


antennas 




(GHz) 


(hours) 


density (m.Jy) 


(mjy/beam) 


50904.92 


1 Apr. 98 


1.28 


VLB A 




1.667 


2.8 


3.6° 


0.23 c 


50908.04 


5 Apr. 98 


4.4 


BR FD HN LA MK NL 


SC Yl 


1.667 


2.1 


221.7 


0.28 


50908.04 


5 Apr. 98 


4.4 


BR FD HN LA MK NL 


SC Yl 


4.987 


2.0 


524.6 


0.59 


50908.04 


5 Apr. 98 


4.4 


BR FD HN LA MK NL 


SC Yl 


15.37 


1.9 


190.6 


0.53 


50978.74 


14 June 98 


75.1 


VLBA +Y1 




1.667 


1.7 


41.5 


0.085 


50985.61 


21 June 98 


82.0 


VLBA +Y1 




1.667 


1.7 


39.3 


0.086 


.Tunc 98 combined 


~79 


VLBA +Y1 




1.667 


3.4 


40.3 


0.075 


50996.65 


2 July 98 


93.0 


VLBA +Y1 




1.667 


2.1 


36.1 


0.084 


51066.45 


10 Sep. 98 


162.8 


VLBA +Y1 




1.667 


6-8 


20.0 


0.036 


51209.44 


31 Jan. 99 


305.8 


MERLIN 




4.994 


10.0 


4.3 


0.15 



a Modified Julian Date (Julian Date - 2400000.5) at mid-point of observations 

b The X-ray flare was first detected on 31.64 March 1998 (MJD 50903.64), and peaked on 1.04 April 1998 (MJD 50904.04). 

c Flux densities for 1 April 98 data after correction for the primary beam, since the VLBA was pointed 11.4 arcminutes from the correct position. 
d 14 June 98 + 21 June 98 data 



NOTE. — VLBA antennas: BR-Brewstcr, FD-Fort Davis, HN-Hancock, KP-Kitt Peak, LA-Los Alamos, MK-Mauna Kea, NL-North Liberty, OV-Owens 
Valley, PT-Pie Town and SC-Saint Croix; Yl - one VLA antenna; MERLIN antennas: Mark 2 (at Jodrcll Bank), Cambridge, Dcfford, Knocking, 
Darnhall, and Tabley. 



were carried out using left and right circular feeds 
with a bandwidth of 32 MHz in each polariza- 
tion centered on 1.7 GHz, except for 5 April 1998, 
which used three bands centered on 1.7, 5.0 and 
15.4 GHz. MERLIN, a 6 clement interferometer 
in Great Britain, observed on 31 January 1999 
at 5.0 GHz, with a bandwidth of 16 MHz in each 
polarization, again using circular feeds. Table 1 
summarizes these observations. The VLBA ob- 
servations were correlated at the VLBA correlator 
in Socorro, NM. Figure 1 shows the epochs of the 
VLBA observations in the context of the radio 
light-curves, taken from the Green Bank Inter- 
ferometer (GBI) (Rupen et al., in prep.; see also 
Clark et al. 2000). The radio spectrum evolves 
smoothly from optically thick (a ~ 0.85, S v cx u a ) 
two days after the X-ray peak, to optically thin 
(a ~ —0.4) synchrotron radiation, 7 days after 
the X-ray peak. CI Cam's radio emission is also, 
for X-ray novae, unusually long-lived, allowing 
it to be observed many months after the X-ray 
flare. The lack of secondary radio flares, in con- 
trast to most other X-ray novae, allows a direct 
link between the radio remnant and the original 
impulsive event. CI Cam is in fact still detected 
(2004) at the ~ 1 mJy level with the VLA (Rupen 
et al., in prep). 

All the VLBA observations were phase refer- 
enced to the strong VLBI sources J0419+5722 
and J0422+5324, while the MERLIN experiment 
was only referenced to J0419+5722. Both of 
these sources are less than 1?5 away from CI Cam 



and are excellent calibrators. For the VLBA 
observations we used a cycle time of 1 minute 
on J0419+5722, 2.5 minutes on CI Cam, then 1 
minute on J0422+5324. The observations were re- 
duced in the standard manner using the NRAO's 
Astronomical Image Processing System (Greisen 
2003). The second epoch was phase and ampli- 
tude self-calibrated, but all other epochs were not. 
During the epoch the total flux smoothly varied 
by 10% in the 1.7 and 5.0 GHz and 20% in the 
15.4 GHz data. The variations were flattened by 
the amplitude self-calibration and did not affect 
the structures in the images. The two June epochs 
were concatenated and imaged together. CI Cam 
evolves slowly enough that this is valid, and it 
improved the image considerably. On 1 April 
the VLBA observed the estimated X-ray position, 
which was off by 11. 4 arcminutes. CI Cam was 
thus at the ~ 60% power point of the VLBA anten- 
nas at 1.7 GHz, so all flux densities for that date 
have been multiplied by 1.7. Since these observa- 
tions were phase referenced the relative positions 
between the high frequency the maps are accu- 
rate, although the relative positions of the 1.7 GHz 
maps may differ by up to a few milliarcseconds 
(mas) due to the ionosphere. Also note that in 
these images lmas= 5d5fc pc AU (d 5 k pc = rf/5kpc, 
with d the distance to the source). A careful polar- 
ization calibration was not done, however we can 
say that CI Cam is not polarized by more than a 
few percent in any epoch, probably less than that 
for the second epoch. 
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Fig. 2.— VLBA 1.7 GHz image of CI Cam on 
1 April 1998, less than a day after the X- 
ray flare. Peak flux density 2.77mJy/beam 
(3.3 x 10 7 K); beam size 9.61 x 3.80 mas at 34?5; 
rms noise 0.23 mJy/beam; contours at ±2 n / 2 x 
0.5mJy/beam, n =0,1,. . . ,4; conversion to bright- 
ness temperature: 12 x 10 6 K/mJy/beam. 

3. Images 

3.1. April 1998 

Figures 2 and 3 show the images from the first 
and second epochs. On 1 April 1998, during the 
radio rise, the source was very faint, and still op- 
tically thick at 1.7 GHz (see Figure 1). The map 
shows a peak of 2mJy/beam; the faint extension 
to the southeast (Figure 2) is probably a noise 
spike. CI Cam may be slightly resolved (4± 2 mas) 
along a position angle of ~ 20°, but this is very 
close to the position angle of the beam (~ 34?5), 
and so may not be real. By 5 April the turnover 
frequency between optically thick and optically 
thin emission was about 5 GHz, and the images 
(Figure 3) show a much more complicated struc- 
ture. At 1.7 GHz, the source is resolved but very 
smooth, extending ~12mas in the north-south 
direction. At optically thin frequencies, 5 and 
15 GHz, CI Cam breaks up into several conden- 
sations with a "kidney-shaped" structure to the 
south. 



3.1.1. Spectral Indices 

There are significant differences between the 
images made at different frequencies on 5 April 
1998: the 1.7 GHz image extends further north 
than the 5 and 15 GHz images would lead one 
to expect, while the western arm of the struc- 
ture seen at 5 GHz is not present in the 15 GHz 
image. Some of these differences could easily be 
instrumental. Interferometers act as spatial fil- 
ters, effectively measuring different components in 
the Fourier transform of the image. The Fourier 
component seen by a single baseline of length B 
corresponds to a spatial frequency of B/X. Thus 
long baselines measure rapid fluctuations of the 
sky brightness, while short baselines respond to 
more gradual changes. The same physical inter- 
ferometer, used at two different wavelengths, may 
respond quite differently to the identical source 
structure. High-resolution interferometric images 
for instance may not recover large-scale structures 
which contain most of the source flux density: one 
sees the crests of the waves, but not the ocean on 
which they rest. For more information see Taylor, 
Carilli, & Perley (1999) or Thompson, Moran, & 
Swenson (2001). 

For CI Cam on 5 April 1998, the 5 GHz VLBA 
image recovers ~ 94% of the integrated flux den- 
sity as measured by the VLA; while at 15 GHz 
the VLBA recovers only ~ 47% of that inte- 
grated flux density. This strongly suggests that 
at least part of the 5 GHz structure which is miss- 
ing from the 15 GHz image has been removed by 
the more extensive spatial filtering operating at 
15 GHz. To test this idea we put the CLEAN 
component model 1 of the 5 GHz data through the 
(u,v) (Fourier) coverage at 15 GHz, to simulate 
what the 5 GHz source would look like if observed 
at 15 GHz. This produces an image with a mor- 
phology almost identical to that actually observed 
at 15 GHz (see Figure 4), and with a total flux 
which is 50% less than the same area in the 5 GHz 
image (i.e., not including the western arm). In 
other words, even if its intrinsic structure were 
identical to that seen in the 5 GHz image, the 
15 GHz image would not show the western arm 
because of the spatial filtering of the 15 GHz ob- 

CLEAN is a deconvolution algorithm which models the 
emission as a set of point sources; see, e.g., Cornwcll 
(Braun). 
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Fig. 3. — VLBA images of CI Cam observed on 5 April 1998. a) 1.7 GHz image: peak flux density 
154mJy/beam (2.4 x 10 9 K); beam size 6.58 x 4.28 mas at 34?5; rms noise 0.28mJy/beam; contours at 
±2 n / 2 x 0.45 mJy/beam, n —0,1, ...,14; conversion to brightness temperature: 16 x 10 6 K/mJy/beam. 
b) 5.0GHz image: peak flux density 131mJy/beam (3.4 x 10 9 K); beam size 1.70 x 1.12 mas at — 14?3; 
rms noise 0.59mJy/beam; contours at ±2™/ 2 x 2mJy/beam, n =0,1,. . . ,12; conversion to brightness tem- 
perature: 26 x 10 6 K/mJy/beam. c) 15.3 GHz image: peak flux density 40mJy/beam (0.35 x 10 9 K); beam 



size 1.02 x 0.58 mas at 18?8; rms noise 0.53mJy/beam; contours at ±2™/ 2 x 2, mJy/beam, n =0,1, 
conversion to brightness temperature: 8.7 x 10 6 K/mJy/beam. 
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servation. In addition, the observed 15 GHz flux 
density would be low by about 50%, in agreement 
with the actual observations. The only notable 
distinction which remains is a ~ 14mJy spur to 
the northeast of the main emission in the 5 GHz 
image, which is not seen at 15 GHz. The spec- 
tral index of this spur must therefore be some- 
what steeper than the a — —0.38 ±0.04 (S v oc v a ) 
which characterizes the rest of the kidney-shaped 
structure. For comparison, the overall 5/15 GHz 
spectral index on this date, as measured with the 
VLA (Rupen et al., in prep.), is —0.36. 

CI Cam is significantly more extended at 1.7 GHz 
than at the higher, optically-thin frequencies on 
5 April 1998. This is evident in Figure 3, but is 
shown more clearly in the super-resolved 1.7 GHz 
image in Figure 5a, where the CLEAN component 
model is restored with a 3 x 3 mas beam. Follow- 
ing the same procedure as above, Figure 5b shows 
the 1.7 GHz model put though the 5 GHz (u,v)- 
coverage. While the southern, kidney-shaped 
emission is fairly similar, the northern component 
is seen only at 1.7 GHz, implying a very differ- 
ent spectral index. Quantitatively, the spectral 
index of the southern, kidney shaped component 
between 1.7 and 5 GHz is a = +0.78 ± 0.03, while 
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Fig. 4. — Contours: the CLEAN model of 
the 5 GHz data, put through the (it, v) cover- 
age of the 15 GHz observation (peak flux density 
79.0mJy/beam; beam size 1.18 x 0.73 mas at 19?6; 
contours at ±2™/ 2 x 7mJy/beam, n =0,2,4,6). 
Colors: the observed 15 GHz image (peak flux den- 
sity 47.2 mJy/beam). 
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that of the northern component must be steeper 
than —1.0. This is most readily interpreted as 
an opacity gradient across the source, with the 
northern component being much more optically 
thin. 

3.2. June 1998 January 1999 

Images from further VLBA observations on 14- 
21 June, 2 July, and 10 September 1998 are dis- 
played in Figure 6. Since CI Cam does not evolve 
rapidly the 14 and 21 June data sets were concate- 
nated, significantly improving the resulting image. 
Figure 7 shows a MERLIN image from 31 Jan- 
uary 1999. The crosses in these images indicate 
the position of the peak of the 1.7 GHz image from 
5 April 1998. Several results are obvious from 
these figures. First, the source does not bear any 
resemblance to a relativistic jet. Second, CI Cam 
has significantly expanded, from ~ 12 mas on 5 
April 1998, to at least 350 mas in the north-south 
direction on 31 January 1999, 306 days after out- 
burst. Third, the expansion is roughly symmetric 
about the center. Fourth, there is significant emis- 
sion within the source's outer rim (shell). Finally, 
the shell's appearance changes significantly with 
time, although the position angle of the major axis 
is constant. In early April, the southern part is 
the brightest; by June- July the northern portion 
is brighter, and the structure looks almost like a 
double ring. In September 1998 the source seems 
more nearly elliptical, and an area in the south 
has brightened from 1.9 mjy in June to 2.7 mjy in 
September, measured at the same position angle 
on the source and approximately the same area. 
Unfortunately, the MERLIN resolution is not suf- 
ficient to see any morphological details, but the 
west side is brighter than the east, showing that 
the source remains asymmetric through at least 
31 Jan. 1999. It is also interesting to note that the 
position angle of the slightly resolved component 
in the first epoch (~ 20°) is comparable to that 
of the major axis in these later images, suggesting 
that this position angle was established within a 
day of the outburst. 

4. Discussion 

4.1. Emission Mechanism 

The radio emission in CI Cam must be primar- 
ily synchrotron, from relativistic particles spiral- 



ing around strong magnetic fields. This is based on 
the high brightness temperatures (up to a few bil- 
lion Kelvin) seen in the VLBA images, and on the 
radio spectrum at high frequencies, which follows 
a steeper power law than allowed for thermal emis- 
sion. Further, any thermal plasma hot and dense 
enough to produce the observed radio emission at 
late times would necessarily produce strong X-ray 
emission as well, at levels well above what is actu- 
ally observed (Parmar et al. 2000; Orlandini et al. 
2000). Interestingly, the peak brightness temper- 
ature in these resolved, optically-thick images is 
well below the canonical 10 11 — 10 12 K brightness 
temperature limit thought to be typical of other 
synchrotron sources, such as extragalactic jets 
(e.g., Kcllcrmann & Pauliny-Toth 1969; Read- 
head 1994) and radio supernovae (e.g., Wicringa, 
Kulkarni, & Frail 1999). Taking the peak of the 
5 April 1998 5 GHz image as an example, assuming 
energy equipartition between relativistic electrons 
and the magnetic fields, and taking the ratio of 
proton to electron energies as k = 100 (as in local 
cosmic rays), the magnetic field is l-2rf 5fe p C G, the 
total energy in relativistic particles and fields is 
1.5 x 10 41 dg^,p 7 c erg, and the synchrotron lifetime at 

5 GHz is 72(£?£.p C days. The energy of the electrons 
producing synchrotron radiation at a frequency v 
is then 304kp C (^/ 5 GHz) 3 / 7 MeV. These values 
should be taken as indicative rather than defini- 
tive. Nevertheless it is clear that, as usual, the 
energy in relativistic particles and fields is not an 
important part of the overall energy budget — the 
total X-ray emission in the flare for instance was 
<~ 9.3 x 10 42 dg kpc erg, while the total integrated 
luminosity was a factor ~ 600 higher (Frontera et 
al. 1998). The rapid radio decay further shows 
that synchrotron losses do not dominate the ob- 
served radio light curve, another not-unexpected 
result. 

4.2. Absorption Mechanisms 

At early times the radio remnant is optically 
thick, indicating substantial absorption. The 
rapid changes in optical depth show that this 
absorption is associated with the source itself, 
rather than some unrelated medium. There are 
two main processes which are relevant here: syn- 
chrotron self-absorption (SSA) by the same rela- 
tivistic electrons which produce the emission, and 
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Fig. 5. — a) The 1.7 GHz image, restored with a 3 x 3 mas beam: peak flux density 1.18mJy/bcam; con- 
tours same as in Figure 3a. b) Contours: the CLEAN model of the 1.7 GHz data, put through the (it, v) 
coverage of the 5 GHz observation (peak flux density 78.0mJy/beam; beam size 2.03 x 1.21 mas at — 10?3; 
contours at ±2™/ 2 x 0.5mJy/beam, n =0,2,. . . ,12). Colors: the observed 5 GHz image (peak flux density 
131.0 mJy/beam). 
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Fig. 6. — 1.7 GHz VLBA images of CI Cam observed months after the X-ray flare. The red crosses mark 
the position of the peak of the 5 April 1998 1.7 GHz image, a) Combined dataset of 14 and 21 June 
1998: peak flux density 4.37mJy/beam (8.7 x 10 6 K), beam size 16.13 x 13.64 mas at 34?5, rms noise 
0.075 mJy/beam; contours at ±2™/ 2 x 0.2mJy/beam, n =0,1,. .. ,6; conversion to brightness temperature: 
2.0 x 10 6 K/mJy/beam. Image made with a Gaussian (u, v) taper of 25 MA (at 30%) and natural weighting, 
b) Image from 2 July 1998: peak flux density 1.45mJy/beam (5.6 x 10 6 K), beam size 11.36 x 10.03 mas at 
13?3, rms noise 0.084 mJy/beam; contours at ±2"/ 2 x 0.2mJy/beam, n =0,1,2,3,4; conversion to brightness 
temperature: 3.9 x 10 6 K/mJy/beam. Image made with a Gaussian (u, v) taper of 25 MA (at 30%) and natural 
weighting, c) Image from 10 September 1998: peak flux density 0.734mJy/beam (1.7 x 10 6 K), beam size 
14.49 x 13.29 mas at 45?4, rms noise 0.036 mJy /beam; contours at ±2™/ 2 x O.lmJy/beam, n =0,1,..., 5; 
conversion to brightness temperature: 2.3 x 10 6 K/mJy/beam. Image made with a Gaussian (u,v) taper of 
15 MA (at 30%) and natural weighting. 
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Fig. 7.— 4.99 GHz MERLIN image of CI Cam 
on 31 January 1999. Peak flux density 
2.159 mJy/beam (4.6 x 10 3 K), beam size 204.69 x 
112.70 mas at 0?3, rms noise 0.15mJy/beam; con- 
tours at ±2"/ 2 x OAmJy/beam, n =0,1,2,3,4; 
conversion to brightness temperature: 2.1 x 
10 3 K/mJy/beam. Image made with a Gaussian 
(u, v) taper of 2 MA (at 30%) and natural weight- 
ing. The black cross, as in Figures 2, 3a and 6, 
marks the position of the peak of the 5 April 1998 
1.7 GHz image. The alignment should be good to 
a few milliarcseconds. 

free-free absorption by dense ionized gas along 
the line-of-sight. The equipartition magnetic field 
derived above would produce an SSA cutoff fre- 
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quency of ~ 3.0c? 5 fepc GHz; below this frequency 
the emission should be optically thick, with a 
characteristic v 2 b spectrum. This is consistent 
with the observed 1.7 and 5 GHz flux densities on 
5 April 1998, which imply a turnover frequency 
around 2.4 GHz. The northern component, seen 
only at 1.7 GHz, could easily be optically thin on 
the same date, since its much lower flux density 
yields a predicted SSA turnover below 1 GHz. At 
the same time the emission from the main, south- 
ern region must truly be as smooth as it appears 
in these images: any significant clumping would 



lead to brightness temperatures above ~ 10 10 K, 
corresponding to SSA turnover frequencies above 
5 GHz. 

CI Cam is also known to have a dense stel- 
lar wind, which could lead to substantial free- 
free absorption. To produce unit opacity at 
3 GHz requires an emission measure (J n^dl) of 

- 3.2 x 10 7 (T e /10 4 K) 1 - 35 pccm- 6 , where T e is 
the electron temperature, n e is the electron num- 
ber density, and I is the line-of-sight path length. 
Taking an n e oc r~ 2 wind, assuming cosmic 
abundances, and taking the inner radius of the 
absorbing material to be 10 AU (corresponding 
roughly to the size of the 5 April 1998 emis- 
sion region), this gives a density at 10AU of 

- 3 x 10- 18 (T e /10 4 K) a675 gcm- 3 . For com- 
parison the modeled density for the wind in 
the sgB[e] star HD 87643 at 10 stellar radii is 
10~ 16 — 10" 15 g cm -3 and 10~ 12 g cm" 3 for the hot 
diffuse and cold dense winds, respectively (Oud- 
maijer et al. 1998). Free-free absorption could 
thus be quite significant, depending on the relative 
orientation of the winds and the radio-emitting re- 
gion, and on the temperature and ionization of the 
winds. In sum, synchrotron self- absorption must 
be important in the early radio remnant, but free- 
free absorption may also play a role. A detailed 
analysis of the multi-frequency radio light curve 
may break this degeneracy. 

Even without any absorption, the radio rem- 
nant shows significant spectral variations: the 
northern component seen at 1.7 GHz on 5 April 
1998 must for instance have an intrinsic emission 
spectrum (f -1 or steeper) which is much steeper 
than the rest of the source (i/~ 0A ). Such extreme 
variations have not been seen in other shell-like 
sources, such as radio supernovae (Bietenholz et 
al., priv. comm.) and supernova remnants (e.g., 
Wright et al. 1999; Delaney et al. 2002). The 
explanation may lie in the very complex environ- 
ment around CI Cam, which could plausibly lead 
to very different densities and magnetic fields in 
different parts of the circumstellar medium. 

4.3. Morphology 

The basic appearance of CI Cam's radio rem- 
nant (edge-brightened, clumpy, shell-like) strongly 
suggests a shock plowing into a dense circumstellar 
medium. The radio emission from all other X-ray 
transients has been either unresolved, or consis- 
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Fig. 8. — Three possible models for the three-dimensional structure of CI Cam. (a) An ellipsoid, where 
patches of the surface are lit up; (b) the ends of a two-sided cone (note that these must both move out and 
widen with time); (c) two spheres. Both the cone and the spheres, for some inclinations, could produce the 
appearance of overlapping circles in later images (Figure 6) . The ellipsoid would require some patchiness in 
the surface, analogous to the brightness variations seen around the rim. The observer is to the right. 



tent with a highly collimated one- or two-sided jet; 
but there is certainly no sign of a jet in these im- 
ages, and the continued brightness of the remnant 
despite substantial expansion requires on-going 
particle acceleration. The relative brightness of 
various features changes substantially from epoch 
to epoch, with some parts of the remnant even 
brightening despite the overall decay; but the over- 
all shape remains fairly constant (see below) . The 
basic implication is that we are seeing the evolv- 
ing results of an explosion inside a dense, clumpy 
environment, much like a radio supernova. The 
changing appearance then corresponds to changes 
in the local shock strength, the particle accelera- 
tion efficiency, or the synchrotron emissivity; these 
in turn are presumably functions of the local par- 
ticle density and magnetic field strength and ori- 
entation. 

Even in our last image, the remnant shock 
is only ~ 350 mas (1750 d^kpc AU) across, and 
is thus still running into material dominated by 
CI Cam's stellar wind. The highly asymmetric 
and clumpy radio morphology probably implies a 
similarly anisotropic and inhomogeneous circum- 
stellar medium (CSM), which seems more likely 
than a severely asymmetric explosion or magnetic 
field, and is consistent with other indications of a 
highly anisotropic wind (Robinson, Ivans, & Welsh 
2002; Hynes et al. 2002). The idea of a dense 
surrounding medium is also supported by X-ray 
observations of CI Cam in quiescence, which show 



strong and highly variable absorption (Parmar et 
al. 2000; Boirin et al. 2002), and by CI Cam's 
strong infrared excess (Dong & Hu 1991; Allen 
& Swings 1976). In this picture the changing 
appearance of the remnant reveals the variations 
in the density of the CSM. The initial bright kid- 
ney shaped structure to the south (Figure 3b and 
c) was presumably due to a particularly dense 
clump on that side, which has been over-run by 
the time of the June- July images (Figure 6a and 
b); while the initially fainter emission to the north 
indicates a lack of dense material on that side. 
Robinson, Ivans, & Welsh (2002) and Hynes et al. 
(2002) suggest that the outflow from CI Cam con- 
sists of a slow (~ 30km/s), dense equatorial wind, 
and a fast (<~ 1000 km/s), less dense polar wind. 
The shock seen in the radio emission between 
2 July 1998 and 10 September 1998 then traces out 
250/8200 years of stellar mass loss, during which 
time the density of the stellar wind changed from 
most dense in the north-west to most dense in 
the south-west; while the similarity between the 
June and July maps shows that the wind did not 
vary dramatically on shorter timescales (75/2500 
years). Significant stellar wind variations are not 
unprecedented. Optical images of ICR+10216, for 
instance, reveal discrete incomplete rings of dust, 
which imply that the mass loss rate of the star 
changes on a hundred year timescale, and is not 
spherically symmetric (Mauron & Huggins 1999). 
Despite these changes in brightness across the 
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source, the overall structure is quite similar from 
epoch to epoch: that is, the brightest regions move 
outward at the same rate as the dim ones. This 
suggests that the radio synchrotron emission sim- 
ply serves to 'light up' the edges of the shock, 
while remaining energetically unimportant. By 
itself this is not so surprising, as the luminosi- 
ties are so low, but this further implies that the 
particle acceleration, and whatever causes it to 
vary around the source, is also not dynamically 
important. This is a common result for shocks. 
The radio remnant of SN 1993 J, for instance, ex- 
hibits changing 3:1 brightness asymmetries around 
a shell which decelerates but remains highly circu- 
lar (Bietenholz, Bartel, & Rupcn 2001; Bartel et 
al. 2002; Bietenholz, Bartel, & Rupcn 2003) and 
the supernova remnant SN 1006 has huge bright- 
ness asymmetries around a similarly circular rim 
(Roger et al. 1988). 

4-3.1. CI Cam in three dimensions 

As with most astronomical imaging, these data 
show a two-dimensional projection of a three- 
dimensional object. The radio remnant is clearly 
elongated, so the intrinsic structure must also be 
asymmetric, either due to an off-center explosion, 
or to an asymmetry in the surrounding medium. 
One puzzle is that the source retains a similar 
appearance while slowing down considerably (see 
below). Several possible geometries, which could 
explain the overall appearance, are shown in Fig- 
ure 8. 

An ellipsoid could explain the outer shape of 
the remnant, with the brighter regions in the mid- 
dle caused by some minor asymmetry in the par- 
ticle acceleration/emission process along the line- 
of-sight to the center. However, it is difficult 
to retain the axis ratio of such an ellipsoid over 
time. Any three-dimensional structure confined 
by denser material around the 'waist' would be ex- 
pected to become more and more elongated over 
time, as the material in the waist should deceler- 
ate more rapidly than that at the poles. Edge- 
brightened cone rims could result from a shock 
lighting the densest regions of a dense wind, as 
in some planetary nebulae, while a pair of spheri- 
cal bubbles might arise in an explosion in a highly 
asymmetric environment. Figure 9 shows a tilted 
circular ring and a two-sided tilted cone that were 
fit to the data; note that the rims of the cone 



both move outward and expand with time. In 
both cases the expansion rate is ~ i - 8 , as dis- 
cussed in the next section. The angle between the 
normal to the titled ring and the line-of-sight is 
53° ± 5°, while the angle between the cone axis 
and the line-of-sight is 23° ± 3°. These fits are il- 
lustrative rather than unique. The two-sided cone 
seems to fit marginally better as it takes into ac- 
count the waist of the nebula. The results for a 
pair of spherical bubbles would be similar to those 
for the two-sided cone. 

4.4. Expansion 

The expansion rate of CI Cam can be fit re- 
gardless of the model used to explain the mor- 
phology, assuming that the source remains self- 
similar throughout. Figure 10 shows the results 
of fitting an ellipse to each image; the axis ratio 
is constant at 1.9:1, and this plot shows the ex- 
pansion of the major axis as a function of time. 
The shock is clearly decelerating, presumably be- 
cause of material swept up from the wind as the 
shock expands. The average radial expansion 
rate of the major axis between an assumed ex- 
plosion date of 31.64 March 1998 (corresponding 
to the first X-ray detection) and 5 April 1998 is 
~ 1.4 mas/day = 11, 000d 5 k pc km/s. This means 
during the 5 April 1998 observations the source 
expanded by ~ 0.5 mas. Since the beam is big- 
ger than this for the 1.7 and 5 GHz images this 
would have no effect. However ~ 0.5 mas is com- 
parable to the 15 GHz beam and therefore the mo- 
tion may have smeared the image slightly in the 
direction of motion. This is probably insignifi- 
cant though, since the structures in this image are 
> 3 mas. Between the first X-ray detection and 
the MERLIN image (January 1999) the average 
expansion is ~ 0.6mas/day = 5, OOOd^kpc km/s. 
A power-law gives a reasonable fit to the ex- 
pansion data: 6 = (4.2 ± 0.9) (t - 50904.09 ± 
0.64)°- 77±0 - 04 , where 9 is the major axis in mas, t 
is the Modified Julian Date (MJD), and the error 
bars are la assuming normally-distributed errors 
and forcing xt ~ 1 • The implied radial velocity is 
v = 14,000(i- 50904.09)-°- 23 km/s. While only 
poorly determined, the explosion date is consistent 
with the beginning of the X-ray flare (on MJD 
50903.64), and the fit is consistent with CI Cam 
being marginally resolved in the first epoch (0 = 
3.6 mas). This expansion rate (t 77 ) is much faster 
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Fig. 9. — Images from (a) June 1998, (b) July 1998, (c) September 1998, and (d) January 1999, overlaid by 
models of a tilted ring (cyan) and two-sided cone (orange). The contour levels are the same as in Figures 6 
and 7. The parameters of the tilted ring model are: radius oc t 8 , angle between normal and line-of-sight 
53° ±5°, and position angle 15° east of north. The parameters of the tilted two-sided cone model are: radius 
oc t , angle between cone axis and line-of-sight 23° ± 3°, cone opening angle 31°, and position angle 7° east 
of north. 



in than the classic Sedov phase seen in supernovae 
remnants (t 0A ), but is consistent with young su- 
pernovae (e.g., for SN 1993J R ~ i°- 92 -°- 78 [Bartel 
et al. 2002]). 

Given the density of CI Cam's circumstellar 
medium, the initial impulsive energy injection 
must have been quite strong to produce such a 
fast shock. One quantitative estimate comes from 
the changing X-ray absorption column seen at 



early times (Belloni et al. 1999). This implies 
column densities N of a few x 10 22 cm~ 2 which 
decline by more than a factor of ten over the 
first two days. It seems reasonable to associate 
this decline with the forcible removal of mate- 
rial by the shock seen in the radio images. By 
the time the X-ray absorption has decayed away 
on 2 April (Belloni et al. 1999) the shock was 
~ 5.4 mas = 4 x 10 14 d5fc pc cm across, according 
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to the above fit. The total number of atoms 
swept up was ~ Nr 2 ~ 8 x ^ 51 d 2 kpc , giving 
an average density of ~ N/r ~ 10 8 d^ pc cm -3 , 
or ~ 3 x 10~ 16 d^ pc g/cm 3 , assuming an average 
mass per atom of 1.4 times the mass of a pro- 
ton. These values are well within the range im- 
plied by spectroscopy during the outburst (< 10 5 
to > 10 10 cm- 3 : Clark, Steele, Fender, & Coe 
1999; Orlandini et al. 2000; Robinson, Ivans, & 
Welsh 2002; Hynes et al. 2002; Miroshnichenko, 
Klochkova, Bjorkman, & Panchuk 2002), and not 
inconsistent with models of the winds in sgB[e] 
stars (e.g., Oudmaijer et al. 1998). The total 
mass swept up by the shock was then of order 
2 x 10 28 dg kpc g, giving a corresponding kinetic en- 
ergy (at 13, 500d 5kpc km/s) of - 2 x I0 i6 di kpc erg. 
This number is very uncertain, given all the un- 
derlying assumptions. Even so it is clear that the 
kinetic energy in the shock could easily match or 
exceed the integrated luminosity of the outburst 
across all observed wavelengths (~6x ^ ib d\ kpcl 
Frontera et al. 1998). 




MJD-50903.64 



Fig. 10. — The expansion of the major axis of 
CI Cam versus time, where time is in days after 
the first X-ray detection on MJD 50903.64 . The 
major axis is estimated by fitting ellipses to the 
images; the errors are taken from the quality of 
that fit, scaled down by a factor ~ 6 to give x 2 ~ 1 
for the power-law model. The five early points are 
taken from the 1.7 GHz VLB A images shown in 
Figures 2, 3a, and 6. The latest point is based 
on the 5 GHz MERLIN observation shown in Fig- 
ure 7. The line is a least squares fit to the data, 
9 = (4.2 ± 0.9)(t - 50904.09 ± 0.64)°- 77±0 04 , with 
O in milliarcseconds and t in days. 

The continued slow deceleration of the shock 



can be used to place constraints on the dis- 
tribution of the circumstellar material via the 
shock jump conditions, in a manner analogous 
to (Chevalier 1982) 's treatment of the expansion 
of radio supernovae. Actually what is derived is 
a relationship between the density structure of 
the outer regions of the ejecta, and that of the 
CSM. Here the 'ejecta' should be thought of as 
the material swept up by the time of the first 
observation; while this is the stellar envelope in 
the case of a supernova, for CI Cam it is probably 
more appropriate to think of dense material sur- 
rounding (accreting onto?) the compact object. 
The shock between ejecta with outer density pro- 
file p e j oc R~ n and a CSM with density profile 
pes ex R w iU expand at a rate R oc t m , with 
m — (n — 3)/(n— b). The observed m = 0.77 
is well within the range observed for radio su- 
pernovae (e.g., Bartel et al. 2002), and suggests 
n ~ 6.3 for an R~~ 2 wind profile. More gener- 
ally one can argue that the observed deceleration 
implies that the shock has swept up a significant 
amount of mass, compared to that of the original 
ejecta. The factor two drop in velocity at late 
times implies a factor two to four increase in the 
mass over the same interval, assuming momentum 
or energy conservation, respectively. 

4-4- 1- Optical Doppler shifts 

Data taken before the April 1998 outburst 
showed a stable spectrum dominated by sym- 
metric, single-peaked hydrogen Balmer lines, 
with wings extending to ±250 km/s (Hynes et 
al. 2002). During the outburst both H I and 
He I lines became much broader and asymmet- 
ric, with line wings extending up to 2500 km/s in 
Ha and He I on 3 April 1998, and perhaps out 
to ~ 5000 km/s in the blue wing of Ha (Hynes 
et al. 2002). Robinson, Ivans, & Welsh (2002) 
also see ~ 2500 km/s line wings in H and He 
on 9 April 1998. While these broad wings faded 
fairly quickly (Hynes et al. 2002), an ultravio- 
let spectrum taken two years after the outburst 
showed P Cygni profiles with absorption wings 
extending to ~ —1000 km/s (Robinson, Ivans, & 
Welsh 2002). It seems natural to associate these 
broad lines with the outflow imaged at radio wave- 
lengths. The proper motions from our images cor- 
respond to a radial expansion of the major axis 
at - U,000d 5kpc (t - 50904.09)~°- 23 km/s. The 
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shock velocity was thus ~ 11, OOOc^ kpc km/s on 
3 April, and - 8, 600d 5kp c km/s on 9 April 1998; 
extrapolating to two years after the outburst gives 
an expected velocity of ~ 3, I00d$kpc km/s. The 
corresponding expansion velocities of the minor 
axis of the fitted ellipse are 60% of these values. 

There are many uncertainties in comparing the 
radio proper motions with the optical Dopplcr 
shifts. First, the source is clearly asymmetric, 
both from the radio images and from the optical 
line shapes at early times, which makes the com- 
parison of radial with transverse velocities some- 
what suspect. Second, while the radio emission 
traces the motions of the outer ejecta, the optical 
lines could be sampling very different regions, with 
different kinematics. This is often seen in super- 
novae (e.g., Bietenholz, Bartel, & Rupcn 2001), 
and is clearly the case for CI Cam at later times, 
when no substantial line wings are seen. Given 
the difficulties the rough agreement between opti- 
cal and radio velocities is remarkable, and suggests 
that the distance to the source is within a factor 
two of the nominal 5 kpc we have adopted. 

4.5. CI Cam as a stellar outflow 

There arc some similarities between the ra- 
dio morphology of CI Cam and the radio out- 
flows/nebulae observed in other stars (e.g., AG Pe- 
gasi [Kenny, Taylor, & Seaquist 1991] and 
HM Sagittiae [Eyres et al. 1995; Hack & Paresce, 
F. 1993; Eyres et al. 2000]), as well as optical out- 
flows around planetary nebulae (PNe) (Balick & 
Frank 2002) and the high mass star Eta Carinae 
(Morse et al. 1998). The hourglass and double 
ring structures seen in many PNe is thought to 
be caused by the density structure of the circum- 
stcllar material into which the ejected material 
expands (Balick & Frank 2002); this is very simi- 
lar to our picture of CI Cam's remnant. However, 
the optical emission from PNe has a very different 
physical origin from the radio synchrotron seen 
in CI Cam, as does the radio emission from most 
stellar outflows (thermal free- free radiation). The 
expansion rate in CI Cam is also much larger than 
the 10 km/s to 1000 km/s seen in other stars (e.g., 
Eyres et al. 1995; Kenny, Taylor, & Seaquist 
1991; Palen et al. 2002). 

The overall structure of the CI Cam radio rem- 
nant (a clumpy shell), its initial expansion speed 
(<~ 10, 000 d5fe pc km/s), and its gradual decelera- 



tion are all very similar to what is seen in radio 
supernovae (e.g., SN1993J: Bartel et al. 1994; 
Marcaide et al. 1995; Bartel et al. 2002 ). There 
are however two basic differences. First, the radio 
shell in CI Cam is decidedly asymmetric, while the 
few resolved radio supernovae are basically circu- 
lar. Supernova remnants on the other hand are 
often asymmetric, which (as in CI Cam) is most 
readily interpreted as the result of a shock run- 
ning into a severely anisotropic medium. Second, 
while radio supernovae - and indeed, virtually all 
resolved, optically-thin radio sources - have rather 
steep radio spectra (a ~ —0.7 to —1.0), CI Cam 
retains a much flatter a ~ —0.4 for many months. 
Since the optically-thin spectral index is a direct 
consequence of the energy distribution of the rela- 
tivistic electrons, this requires a very unusual (but 
also quite stable) particle acceleration process. 

4.6. CI Cam as an X-ray transient 

While radio emission is often associated with X- 
ray transients and binaries (e.g., Hjellming & Han 
1995; Fender & Kuulkers 2001), CI Cam is unique 
in showing a slowly-expanding, two-dimensional 
remnant. In all other cases, any resolved radio 
emission is dominated either by separable ejecta, 
or by a collimated jet, often moving at relativistic 
speeds (e.g., Mirabel & Rodriguez 1994; Hjellming 
& Rupen 1995; Fomalont, Geldzahler, & Brad- 
shaw 2001; Mioduszewski et al. 2001; Blundell 
et al. 2001; Gallo et al. 2004). This has led to 
the wide-spread assumption that radio emission 
from X-ray binaries always indicates the presence 
of a relativistic jet, and many models now take ra- 
dio emission as the best unambiguous tracer of jet 
emission (e.g., Fender 2003). CI Cam throws this 
entire paradigm into doubt. Does this have wider 
implications for X-ray binaries in general? or is 
CI Cam truly unique? 

4-6.1. CI Cam as a collaspar? 

One possibility is that CI Cam did indeed pro- 
duce a fast jet, but this jet was quickly smothered 
by the dense surrounding material. In this model 
the energy in the jet was converted into a less di- 
rected outflow, which we see as a shock moving 
out into the CSM. This is similar to the failed-jet 
model for supernova explosions in collapsars (e.g., 
MacFadyen, Woosley, & Heger 2001). Quantita- 
tively, one expects a relativistic jet to slow catas- 
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trophically when it has swept up an equivalent rest 
mass energy (e.g., Rhoads 1997; Sari, Piran, & 
Halpcrn 1999): 

E jet = f 7TC 2 (p) r 3 stop 

where Ej et is the energy in the jet, and (p) is the 
average density within the radius r stop at which 
the jet slows catastrophically. This equation is ap- 
propriate to a jet which is stopped quickly enough 
that lateral expansion can rapidly isotropize the 
energy. For CI Cam, we know that the remnant is 
far from jet-like by a radius of 6 mas (30d^ k pc AU) , 
and according to the arguments in §4.4 the energy 
in the shock (in this model derived entirely from 
the original jets) is of order ~ 2 x 10 46 c?5 fcpc erg. 
This gives a lower limit to the mean density re- 
quired to stop the jets of (p) > 10~ 19 c?5 kpc g/cm 3 , 
well within the plausible range for the wind of a 
sgB[c] star. It is interesting to note that simula- 
tions of "stalled" jets in the context of 7-ray bursts 
predict morphologies very similar to what we see 
in CI Cam (see particularly the 36°-inclination im- 
ages in Figures 8 and 9 of Ayal & Piran 2001). 

The question arises as to why other jets are not 
similarly smothered. Most other X-ray transients 
live in relatively pristine environments, where the 
jet might more easily escape. In binaries with 
substantial circumstellar material the accretion 
tends to be more continuous (e.g., Cyg X— 3 and 
SS433), and the resulting jets are also 'on' more 
often, and are thus perhaps capable of sweeping 
away the CSM. This fits in with Hynes et al. 
(2002)'s picture of the CI Cam outburst, in which 
a compact object in a highly eccentric long-period 
(> 100 years) orbit moves through the dense equa- 
torial wind of its companion star; the brief period 
of super-critical accretion, when the compact ob- 
ject is close to the companion, leads to the X-ray 
flare and the expanding radio shock front. A long- 
period orbit is required to explain the lack of pre- 
vious known flares, which would have been obvious 
in the optical to historical astronomers. 

4-6.2. Shock signatures and implications for 
other X-ray transients 

Given that CI Cam's radio remnant is clearly 
dominated by a shock, it is natural to ask whether 
the emission at other wavelengths could also be 
associated with that shock, rather than with the 



accretion disk. Extrapolation of the radio syn- 
chrotron spectrum to other wavelengths shows 
that the expected synchrotron emission is far be- 
low that observed during the outburst (Bclloni et 
al. 1999; Hynes et al. 2002). Similarly the thermal 
emission produced by the shock is likely to be quite 
modest, for three reasons. First, the shock has not 
had enough time to become radiative, for any rea- 
sonable densities. Second, the characteristic tem- 
perature of the forward shock (~ 300 keV for a 
10,000 km/s shock) (a) lies outside the passband 
of most available instruments, and (b) puts the 
shocked gas in a regime where cooling is quite slow. 
Third, while the reverse shock should be slower 
and hence lead to lower and more easily observable 
temperatures, the resulting fluxes arc likely to be 
very low; scaling from the radio/X-ray ratios in ob- 
served supernovae for instance (e.g., Pooley et al. 
2002) the expected emission is orders of magnitude 
below the observed flux during the outburst. The 
observed soft X-rays, in this as in other X-ray tran- 
sients, probably arise in and around the accretion 
disk rather than in the shock. Thus the only obvi- 
ous signatures of the shock outside the radio rem- 
nant are the optical and UV spectral lines. Some 
of these are Doppler broadened by a thousand or 
more km/s; others evolve very rapidly, requiring 
very high densities, as might result from a shock 
(Hynes et al. 2002). Recent spectra (2001-2002) 
still show many spectral lines at very different 
strengths than their very stable pre-outburst levels 
(Hynes et al. 2002; Miroshnichenko, Klochkova, 
Bjorkman, & Panchuk 2002), which suggests some 
continued disturbance around the binary. None of 
these is really a 'smoking gun,' indicating that, for 
an unresolved X-ray transient, a shock and a freely 
expanding jet may not be distinguishable. 

It does not however seem likely that the radio 
emission from many other unresolved X-ray bina- 
ries can be fully explained by shocks similar to 
that seen here. CI Cam's radio light curve is un- 
usual in showing a single outburst followed by a 
smooth, steep-spectrum (S v oc v~® A ) decay. Out- 
bursts from other sources often lead to multiple 
radio events, suggesting strong on- going activity 
more readily associated with the central source 
than with a single expanding shock wave. There 
could be a series of shocks of course, but that re- 
quires re-furbishing the surrounding medium on 
rather fast timescales (days to weeks). On the 



14 



other hand, this single shock evacuating the sur- 
rounding medium may explain why there was only 
one explosive event, i.e., the shock got rid of any 
additional material that could have been accreted 
by the compact object, and the compact object 
moved away from the high density region before 
the stellar wind could replenish the material. The 
quasi-steady radio emission seen in the low/hard 
state is also very different from CI Cam's optically- 
thin remnant, and requires a very small and con- 
tinuously replenished radio source. Given these 
indications of very different power sources, it is 
quite surprising that CI Cam lies nicely on the 
loose X-ray/radio correlation found for X-ray bi- 
naries (Fender & Kuulkers 2001; Gallo, Fender, & 
Poolcy 2003), with a scatter typical of other X-ray 
transients. It is far from obvious why direct syn- 
chrotron emission from jets should have the same 
efficiency as particle acceleration in an expanding 
shock. 

5. Conclusions 

We have presented evidence that the radio rem- 
nant of CI Cam's 1998 outburst takes the form of 
a clumpy, decelerating shell, strongly suggesting 
an expanding shock front, rather than the rela- 
tivistic jets seen in many other X-ray transients. 
We suggest that the reason for this distinction 
lies in CI Cam's very dense stellar wind, which 
could easily convert even quite energetic jets to 
a quasi-spherical outflow. The very broad spec- 
tral lines seen in optical emission early on, and in 
the ultraviolet several months later, seem to be the 
only other unambiguous signatures of this outflow. 
Associating the radio outflow with those Doppler 
shifts suggests that the adopted distance of 5 kpc 
is probably good to a factor of two. Assuming 
the changing X-ray absorption at early times is 
due to circumstellar material removed from the 
line-of-sight by the radio shock, we estimate the 
shock energy as ~ 2 x 10 46 dg fepc erg, comparable 
to the total integrated luminosity of the outburst. 
If the smothered jet model is correct, this repre- 
sents one of the very few measurements of a jet's 
total kinetic energy. The relativistic jets seen in 
other X-ray binaries and transients must either 
be much more powerful, or be expanding through 
much less dense material. The radio light curves of 
most unresolved sources suggest continued emis- 
sion from a small region, so it also seems unlikely 



that extended shocks dominate the radio emission 
for these sources. CI Cam seems to be related to 
other X-ray transients, as collapsar supernovae are 
related to 7-ray bursts. 
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